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Abstract 

We introduce an approach to expand gauge-invariant Wilson operators on 
lattice. This approach is based on non-abelian Stokes theorem and overcomes 
some shortage of some former methods. It is also suitable for expanding any 
Wilson operator on lattice. 
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1 Introduction 



Expanding gauge-invariant Wilson operators, or, Wilson loops, according to lat- 
tice spacing a is a very important work when we do theoretical analysis in lattice 
QCD. It is a necessary step in many cases, such as the discussion of improved 
actionB, the exploration of improved cooling0 or the investigate for renormaliza- 
tion procedure^ when study instanton phenomenon, etc. Furthermore, in order 
to calculate masses of glueball states, we developed a new opinion@, in which the 
perturbative expansion of Wilson loops is also a needed progress. 

Due to the complicated of path ordering operators, it is cockamamie and difficult 
to expand the Wilson loops according to normal method if the number of links of 
Wilson loop is large or it is needed to expand Wilson loop up to higher precision. 
There is a prevalent approach introduced by Luscher and Weisz in ref . . But this 
approach has some shortcoming, for instance, one is likely to face some difficult 
when he tries to expand relatively complicated spatial loops. 

We introduce a practical easy approach which is based on the non-abelian Stokes 
theorem here. This approach has overcome the shortcoming in ref. and is 
therefore suitable for any Wilson loops on lattice. It is also suitable for any of above 
cases in ref. ^ |^, 3. To do this, we will give a short comment on Luscher- Weisz 
method in the following section and some mathematical preparation in section 3. 
We exemplify our method in section 4. 

Since the case in ref. is the most difficult and universal one, we should pay 
our most attention on this case. 

2 The Comment on Luscher- Weisz Method 

Let us observe an arbitrary gauge-invariant Wilson loop on the cubic lattice. In 
order to get the perturbative expansion of the loop up to O(a^) ( To simplify, 
this paper only discusses the expansion up to 0{a^ j), Luscher and Weisz adopt 
the following method: they fixed a gauge at first, then try to rewrite it in gauge- 
invariant form. 

Assume that one needs to expand the plaquette operator: 

= HHP^An) = T.T. Tr{U{n, n)U{n + /i, v)U-\n + P, ^^)U-\n, //)}, (1) 

where fi is the unit vector along ^ direction, U {n, /i) is link variable defined as 
follows: 

1 

U{n, fi) = T exp a J dtA^{an + ap,t) , 


in which A^is gauge potential^. We point out here that the summation is different 
following different problem. For instance, the summation of n is over all the 4- 

^ We use a notation here that Afj_ is anti-herniitian, traceless N x A^-matrix 
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dimensional lattice sites and the summation of directions( fiv) is over all possible 4- 
dimensional space-time directions when one discusses the improved actionE , while 
the summation of n is over spatial lattice sites in the same time slice and that 
of directions is variant in 3-dimensional spatial directions according to different 
problems in the case of ref. ||^. 

One can always choose a proper local gauge on lattice to make the gauge field 
satisfy 

Af^{x) = 0, for all x, 

Au{x) = 0, for all x with = 0. (2) 

pl^lu{0) then reduces to 

1 

pl^u{0) = TrT expa J dtA^{afl + ajlt) . (3) 



At the same time, at x = 0, one has 

a^a^+M, = D^DIF,,. (4) 
Insert Eq. into Eq. (^) and expand Eq. (^) according to lattice spacing a: 

+ }-a'Tr{Fl, + F,,{Dl+^-D,D, + Dl)F,, 

+ -^{D^ + D,)F,,{D^ + D,)F^,} + 0{a'), (5) 

where F^^ = d^A^ - d^A^ + [A^, A,,] and D^F^p = d^F^p + [Ap,F^p]. Due to 
periodic boundary condition, one can ignore all the (covariant) differential terms 
in the last result ( Section 4 shows more details) and gets: 

plpM = \a'TrFl, - ^^a'Tr[Fp,{Dl + Dl)FpJ\ + \a'TTFl + 0{a'). (6) 

Now let us observe Eq. (H). In this case, because 

DpD,-D,Dp = [Fp,,-l (7) 

the corresponding in Eq. @ is unambiguous. But, when one tries to expand 
some complicated spatial Wilson operators, l.h.s. of Eq. is possible as the 
form d^duAp{fi ^ v ^ p). There is ambiguity that such terms may correspond 
both to DpFyp and to DyF^p^ which means one can not determine the expansion 
unambiguously, especially when the summation of directions is not over all possible 
directions but some definite onesS. 
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In the case of ref. Q that one should sum over all possible 4-dimensional 
directions, Luscher and Weisz eliminate the ambiguity through the presumption 
that the expansion of any Wilson loop after summation over all possible directions 
and lattice sites is invariant under any Lorentz transformation( In 4-dimensional 
Euclidean space, it is 5*0(4) group), parity and charge conjugacy transformation. 
This presumption will some time complicate our expansion, furthermore, it has not 
been proved strictly. We show a debate here. Suppose we accepts this presumption 
and restricts ourselves into 3-dimensional space, the subspace of the 4-dimensional 
Euclidean space and time. One natural deduction is that the operator, which has 
been summed over all space lattice sites in the same time slice and all possible 3- 
dimensional directions, is invariant under S0{3)^'^' group since the S0{3) group is 
the subgroup of the 5*0 (4) group. As such operator belongs to representation Af^, 
the trivial representation of the cubic group on lattice, one may regard that the 
J^*-^( in continuum limit) of the operator belonging to Af^ is O"*""*". In other words, 
the subduced representation of the representation 0++ of the S0{3)^'-^ group can be 
reduced into Af'^l subduced representation: the representation which is obtained 
by trivially embedding the cubic group O^*" into the group 5'0(3)^*"). But, as we 
know, there also exist other representations in the group 5'(9(3)^*", the subduced 
representation of which can be reduced into Af^, for instance, 4"*"+ 00, which 
means that the J^'"( in the continuum limit) of the operator belonging to Ai~^ 
is not always So, the validity of the presumption is doubtful. In fact, a 



discussion in ref. [|T0| does not support this presumption too. To say the least, 
even this presumption is right, one can hardly apply it to expand the operator, for 
instance, the PC of which is not ++. 

To solve the problem, we developed a new approach which can determine the 
perturbative expansion unambiguously. The next section is the proper mathemat- 
ical preparation. 



3 Mathematical Preliminary 

For operators, or loops, the traces over color space are not performed in this section, 
so, these operators are gauge covariant under gauge transformation. 

Now we begin our preparation with some simple examples. 

Observe the operator W{0, Z) shown in Fig. 1: 

W{0, Z) = Ur{Z, O; Cx)W{Z)Ur{0, Z- Cx), (8) 

where is a definite gauge covariant operator which is from point Z and back to Z, 
and Ur{Z, O, Cx) = -Pexp dtAi{it) is connector. Some properties of connectors 
are list as follows®: 

a) 

[/h(X3,Xi;Ci + C2) = [/«(X2,Xi;Ci) • [/^(Xg, X2; C2), (9) 

where Ci is a curve going from point Xi to point X2 and C2 is a curve going 
from point X2 to point X3. 
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b) 

UR{X,Y;C)Uii{Y,X;C) = l, 
UUY,X,C) = Un{X,Y;C), (10) 

C is a curve from X to y and C the same curve but oriented in inverse 
direction, running from Y to X. 

From the opinion of differential geometry, W is the shifted version of operator W 
from point Z to point O. 

Under local gauge transform V we have: 

W{0,Z) V{0)W{0,Z)V-\0), 

W{Z) V{Z)W{Z)V~\Z). (11) 

Let W{Z) = W + xdiW-\ l-^d^W-l , where W = W{0) is the same operator 

as W{Z) but it is from O and back to O. 
We can prove that 

W{0, Z) = W + xDiW + ■ ■ ■ + ^D'lW + • ■ ■ . (12) 

Then we study the operator described in Fig. 2 ( We also denote it by W)\ 

W{0, Z) = Ur{Z, O- Cy + Cx)W{Z)Ur{0, Z- Cx + Cy). (13) 

We let W{Z) = W + xdtW + ydiW + \{x'^df + 2xydidjW + y'^d'^^W) + ■ ■ -, in which 
W = W{0). Using Eq. (^) and the properties of connector we get 

W{0, Z) = W + xDiW + yDjW + ^(x^Df + 2xyDjD, + y''D])W + ■■■. (14) 

2 ■' 

It should be pointed out here that, in general, the order of DjDi can not be reversed 
in above equation. 

Other more complicated case can be calculated as similar as above. 

We write out non-abelian Stokes theorem as( We let to = o)0 

U^,,{x) = U^{x)U,{x + Ji)Ul{x + d)Ul{x) 

= Fexp I ds I dtFnu{x + sfi + v) 
Jo Jo 

<x> rn „a r-U-i 
1 + XI n / '^■^W dtiF^^{x + Sn'jl + trfi)---F^^{x + Sip. + 1111). 

(15) 



Here F^y{x + sp + tV) is the shifted version of F^y{x + sfi + W) from point x + sfi + tv 
to point X, or: 

Vis.t) = Fexp / dt'AJx + t'V) - Fexp / ds' AJx + W + s'p.), 
Jo Jo 
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F^^{x + SnfL + tnV) = V{s,t)F^^{x + sJl + tV)V^{s,t). (16) 

Now using Eq. (p^), (|15|), ( |T^ ) we can get the expansion of the operator from O 
back to O which is shown in fig. 3 (a): 

U^j\n) = Ui{n)Uj{n + i)UKn + j)U^^{n) 

= 1 + a'F,, + j{D, + D^)F,^ + j{D^ + D] + \d^D,)F,^ 

+P^, + '^Ff^ + 0{a'), (17) 

where P^- = ^Ff^ + ^Fij(Df + D'j)Fij, Fij = Fij{n). Because our precision in last 
result is O(a^), we have ignored all the full (covariant) differential terms the ranks 
of which are large than or equal to 0{aP){ More details are shown in section 4), for 
instance the term ^Fij{Di + Dj)Fij. 

Using Eq. (|T7D and the properties of connector we get the perturbative expan- 
sion of the operator in Fig. 3(b): 

Ulf{n) = Uj{n)Ul{n-l + ])U}{n~l)U,{n~l) 

= 1 + a'F,, - ^(A - D,)F,, + ^{D^ + D| - ^D^D,)F,, 

+P,, + ^4 + 0(a^). (18) 

With the same reason we get the perturbative expansion of the operator in Fig. 
3(c): 

U^fin) = UKn-t)U]{n-l-m{n-i-j)U,{n-]) 

= 1 + a'F,, - y (A + D,)F,, + ^(A' + + ^A-A)^^.- 

+P,, + ^4 + 0(a^). (19) 
And the expansion of the operator in Fig. 3(d): 

Ug\n) = U^^(n-m{n-j)U,{n + i-])Ulin) 

= 1 + a^F,, + y (A - D,)F,, + j{D^ + - ^AA)^^.' 

+P,, + ^i^5 + 0(a^). (20) 

One can prove that the expansion of the charge conjugacy parity of above operator 
f/(*) is equivalent to replace Fij by -Fij in Eq. (0) - (|2|). 

Now we can expand arbitrary Wilson loop on lattice using above equations. It 
is shown in the following section. 
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4 The Expansion of Arbitrary Wilson Loop on 
Lattice 



We now give a general properties of the perturbative expansion of Wilson loops on 
lattice before we perform the expansion. 

a Because differential operator can be presented by combination of difference 
operators with various ranks, the final contributions of full differential terms 
are zero due to the periodic boundary condition and we shall ignore such 
full differential terms in the expansion. In this reason, one can also expand 
Wilson loops on any lattice site with the same result. 

In this reason, one can ignore full (covariant) terms with rank large than or 
equal to in section 2 and 3 in our case. 

b On lattice an arbitrary Wilson loop can be expanded as 

oo 

0{n) = Y.a^Ok{n), (21) 

fc=0 

where Ok{n), with mass-dimension k, is gauge- invariant polynomial of 
and Ffj,i, and their derivatives. In above equation even k corresponds to 
parity P = + while odd k correspond to parity P = —. 

c As we know, the C parity of the real part of the operator is positive and 
that of the imaginary part is negative. Since the transformation of covariant 
differential Di{Di- = di ■ — •]) can be regard to be l""*" under S'C'(3)^'^ 
group in this case, we can regard Di as a real operator in the expansion. 

One can simplify the expansion using three above properties. For instance, from 
above we know that Pij corresponds to PC = ++ and it is invariant under trans- 
lation on lattice if we only expand operators up to precision O(a^). 

Utilizing the perturbative expansion of U^'\i = 1, 2, 3, 4) in (|13) - one 
can expand an arbitrary Wilson loop on lattice, we exemplify it here. 

The expansion of operators with PC = ++ is the most complicated and we 
should only focus on it. So, when we say an gauge-invariant operator L here we 
actually means ^{1 + C + P + C P) L, where P and C are parity operator and charge 
conjugacy parity operator which acts on operator L. 

Let us observe the operator described in Fig. 4 and only care terms with 
PC = ++, so that we may only keep real terms with even k in Eq. (|2l|) . 

L = Tr(l-2-3-4-5-6-7-8), (22) 

where, for simplification, we denote, for instance, the link variable U{1) by 1 and 
f/t(l) by 1, etc. 
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We rewrite L as 

L = rr[l-2-3-4-9-l-l-9-5-f0-10-6-7-8] 

= Tr[ph-pk-pm, (23) 

where p/i = 1 ■ 2 • 3 ■ 4 ■ 9 ■ 1, = 1 ■ 9 ■ 5 ■ 10 and p/3 = 10 • 6 • 7 • 8. Each of the pi 
is gauge covariant respect to point O under local gauge transformation. 
Using equations in section 3, we get 

3 

ph = l + a''Fik + ^{Di + Dk + 2Dj)Fik 

3 

+ -(^2 ^ ^2 ^ 3^2 ^ ^ 3^.. ^ 3^.^)i7^.^ + p.^ + o(«7)^ ^24) 
where Dik — DiDk, etc. 

p/2 = 1 + a^Fjk + ^(D,- + + j{D] + Dl + ^D,k)Fjk + 0{a'). (25) 

and 

pk = l + a^F,k - y (A - Dk)F,k + ^(D^ + - ^D,k)F,k + O(a^). (26) 
So, up to 0{aP), we have: 

L = 1 + 2a''Tr{FikFik) + -Tr{FjkFjk) + 2a''Tr{FikFjk) 

+^Tr[F,,(4D2 + Dl + 3D] + 6A,)F,,] + ^Tr[F,k{D] + Dl)F,k] 

+^rr[F,fe(2D2 + D] + Dl + 3A,-)^ife] + 0{a'). (27) 

Here we have ignored terms with PC 7^ ++ and full covariant differential terms. 
Of course, we can do as similar for terms with PC ^ ++ in the expansion. 

5 Conclusion 

We show an approach in the paper for the expansion of the Wilson loops on lattice. 
The approach presented here can determine the perturbativc expansion of Wilson 
loop on lattice unambiguously according to P and C and is suitable for any Wilson 
loop on lattice despite how tanglesome the Wilson loop is. 

It is also an appropriate method to discuss the improved action or improved 
cooling on lattice, etc. In fact, in some of the former papers which discuss improved 
action, there maybe exist some terms in higher precision which arc not invariant 
under the rotation in 5*0(4) group as the continuum theory requires, and therefore 
should bring up some artifact on lattice. Using this approach may help us to elim- 
inate those unwanted terms if we discuss the improved action or other interesting 
problem up to higher precision. 
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Figure captions 

• Figure 1 The operator Ur{Z, O; Cx)W{Z)Ur{0, Z; Cx), where Ur is the 
connector described in the paper. 

• Figure 2 The operator Ur{Z, O; Cy + Cx)W{Z)Ur{0, Z; Cx + Cy), where 
the connector Ur{Z, O; Cy + Cx) = Ur{Y, O; Cy) ■ Ur{Z, Y; Cx). 

• Figure 3 The operators C/*^*-* which are from point O and back to point O. 
They are of the same shape with the different starting and ending point. 

• Figure 4 The operator L — Tr{l •2-3-4-5-6-7-8). Here we denote the 
hnk variable U{i){i = 1, 2, • • • , 10) by i and U'^{i) by i. 
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